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We study the chiral Gross-Neveu model with Wilson fermions. In the framework of the 
Schrodinger functional we show that in general not only the bare mass has to be tuned to achieve 
chiral symmetry in the continuum, but also coupling constants multiplying chirally non-invariant 
four fermion terms. The strategy for fixing the parameters of the theory is explained in pertur- 
bation theory and the results of a first order calculation are presented. The results are shown to 
agree with the infinite flavour limit. 
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1. Introduction 

Theories of self-coupled fermions in two dimensions may show interesting features like asymp- 
totic freedom and chiral symmetry [p. A model with these properties provides a lab for all sorts 
of questions raised in the context of lattice formulations of chiral fermions in four dimensions. In 
principle, the universality of different fermion discretizations can be investigated, new ideas can be 
tested in a controlled environment or "established" concepts can be challenged with high precision 
Monte Carlo (MC) data. 

Many of the perturbative and nonperturbative properties of the continuum chiral Gross-Neveu 
model (CGN) have been calculated, see e.g. [§] and I^J. The properties of chiral symmetry with 
Wilson fermions have been studied in the limit of an infinite number of flavours [@, §]. In this work 
we use lattice perturbation theory (PT) to compute the critical mass and other parameters of the 
lattice action with Wilson fermions. 



2. Self-coupled fermions 

Rather than just stating an action associated with CGN we start with a general action contain- 
ing all possible terms that self-couple the fermions. Their number will be reduced by symmetries 
and interdependence. In this way we are sure to include all interaction terms that might be impor- 
tant for renormalisation. 

In two dimensions fermion fields 1/7, y/ have mass dimension 1/2. Therefore the couplings 
multiplying four fermion interactions (4FI) are dimensionless and relevant. In the continuum eu- 
clidean action for N flavours of Dirac fermions 1 

5 = / ^ {^"^ ((T^)«v3^ax -f" V^3, -I- V^«/ V^sy V^5/ } (2.1) 

the 4FI consist of four Dirac spinors with Dirac (Greek letters a,j8, . . .) and flavour (Latin letters 
/, j, . . . ) indices contracted by the tensor C. 

Imposing now a set of symmetries: flavour U(N) symmetry, Euclidean invariance, parity and 
using the properties of /-matrices in two dimensions (7^/5 = i^vYv) one finds that C can be written 
as a superposition of 6 terms characterized by three different Dirac structures and by two different 
flavour structures: 2 

Oss= (W) 2 o' ss ^ a m a wf 

Opp = OFow) 2 o' PP = l«(W5A>) 2 (2.2) 

Ow = vO 2 o' vv = L^ a (wr^ a v) 2 

In the flavour non-singlet terms on the right hand side the sum is over the N 2 — 1 generators X a of 
SU(iV). But these terms are not all independent. Fierz transformations (FT) allow to interchange 



products of Dirac spinor bilinears like the terms in ( |2.2| ). In two dimensions the FT yield three 
equations like 

OsS — Opp —jjOyy — \0' vv . (2.3) 



x If not stated otherwise, repeated small letter indices are summed over. 

2 The contraction of Dirac and flavour indices is obvious, hence they are suppressed. 
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Taking these interdependencies into account there are three independent 4FI and hence three di- 
mensionless couplings 3 . 

Chiral symmetry is not considered at this stage, because we want to study the theory on the 
lattice with the Wilson discretization, which breaks chiral symmetry explicitly. One can choose 



any three terms from ( |2.2| ). In view of nonperturbative MC calculations we choose the three terms 
with flavour singlet bilinears only. To be specific, on a cubic lattice with extension T x L we study 
the theory with the euclidean action 

S = al Z \W(D W + m ) y - X - (g 2 s O ss + g\ PP + g\ O vv ) j , (2.4) 

where D w is the Wilson-Dirac operator and a the lattice spacing. 

3. Chiral Ward-Takahashi identity 

The PC AC (partially conserved axial current) relation of the theory obtained in the naive con- 
tinuum limit of ( pT4t > takes the usual form, but is modified by a term originating from the scalar and 
pseudo scalar 4FI 

(0(y)d^(x)) = (0(y)2 m P(x))-(0(y)(g 2 s -g 2 P )/2P(x)S(x)), x^y, (3.1) 

where 

Ap(x) = W(x)YnT^¥(x) , P(x) = \j/(x)y 5 \i/(x) , S(x) = . (3.2) 

The vector 4FI is chirally invariant on its own. Note that an U(l) chiral transformation was used 
to derive this and therefore only the flavour singlet current and densities appear. Imposing chiral 
symmetry in the continuum theory the mass has to vanish mo = and the couplings of the scalar and 
the pseudo scalar 4FI have to be equal gp = g 2 s . The PC AC relation then turns into the conservation 
of the axial current. The resulting action with the interaction — \ (g| (Oss + Opp) + gy Ovv) is 
what in the literature is called the CGN 4 . 

However, at finite lattice spacing chiral symmetry is explicitly broken by the Wilson term, the 
mass is not protected against additive renormalisation and gj, = g 2 s is no longer protected. The 
right hand side of ( fTI! ) becomes a mess due to operator mixing. To achieve chiral symmetry in the 
continuum limit means for properly normalized fields and axial current 

(0(y)d tl (A tl ) R (x))=0(a), (3.3) 

where is a symmetric lattice difference operator. That is in an asymptotic expansion of correla- 
tion functions (0(y) (A /i ) R (jc)) in the lattice spacing the zeroth order vanishes. This condition 



is used to fix the parameters of the action ( j2.4| ) in MC simulations [|6p and in PT. 

4. Correlation functions 

We define correlation functions in the Schodinger functional (SF) setup ^ in order to uti- 



lize ( p.3| ) [|8p. In this finite size scheme the fermion fields obey Dirichlet boundary conditions in 



3 For N = 1 the number of degrees of freedom allow only for one local interaction term. Here we consider N > 1. 
4 In most large N studies the vector term is not present. 
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Figure 1: Tree level diagrams for fx. 



time and generalized periodic boundary conditions with a phase shift parametrized by an angle 6 
in space 

\jf(x + L) = e'VO) , W(x + L) = Q- w \j/(x) , e [0,27r) . (4.1) 

The Dirichlet boundary conditions in time may give rise to new counter terms defined on the bound- 
ary, in the sense that the associated bare coefficients are needed to absorb infinites on the way to 
the continuum limit. Relevant operators or composite fields living at the boundary are those which 
have mass dimension one or less. Thus objects like y/Ty/ with r = {1, /o, /i, 75} can appear. But 
all these terms either vanish because they violate parity or for the boundary fields that we employ. 
Thus no additional boundary counter terms have to be considered. 
The correlation functions 

fx(xo) = -^- £ (xj/(x)r x¥ (x)X( yi )Y 5 azi)), r A = 7 o75, r> = 75 , (4.2) 

to be considered are correlators of a zero momentum pseudo scalar boundary state built from the 
boundary fields £, £ and insertions of the time component of the axial current (/a) and the pseudo 
scalar density (//>) respectively. For details on the definition of the generating functional and its 
perturbative expansion we refer the reader to a forthcoming publication 



5. First order renormalization 



In order to calculate the parameters that ensure (3.3) in first order of PT, the correlation func 



tion (4.2) has to be expanded to this order 



fx(xo) = A°\xo)+g 2 s f^ S \xo)+g 2 pf^ P \xo) +g 2 vf^ V) (x )+O(g 4 ) , (5.1) 

where we assume all three couplings to be small and comparable to some small g 2 , i.e. 

gj = c ig 2 , and Q = 0(1), i = S,P,V. (5.2) 

The tree level amplitude f^ is the sum of the two diagrams sketched in Fig. [l[ In these diagrams 
the dotted lines represent the time slices xo = and xo = T. Plain lines are used for the fermion 
propagator and the small open circle in the middle symbolises the insertion. 

The first order amplitude f^ = csfx^ +cpfx' P \xo) +cyfx^ is a sum of three contri- 
butions from the three 4FI. Each such contribution is itself a sum of 10 diagrams. Due to the 
75-hermiticity of the Wilson Dirac operator, and thus the fermion propagator, the diagrams la and 
lb and also 2a and 2b in Fig. |2]are equal. In these diagrams the small filled circles (dots) represent 
the 4FI. The diagrams in the lower row of Fig. || are very similar to the respective ones in the upper 
row but have two separated dots. Because the 4FI consists of two Dirac spinor bilinears, there are 
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Figure 2: First order diagrams for fx. 



always two similar but distinct contractions: either the flavour and Dirac indices of the four spinors 
are contracted collectively or the two bilinears are contracted separately. One single dot refers to 
the former and two separated dots to the latter possibility. 

We define now the renormalised correlation function (f A (xo)) R and h A (0,xo/L) as dimension- 
less quantities on a lattice with fixed ratio T/L 



(/aOo))r = Z A Z% f A (xo) and h A (e,x /L) = Ld (/aOo))r , 
where we introduced normalization factors for the axial current and the boundary fields 



z A = 1 + izf + 0(g 4 ) , z c = 1 + g'Z^> + 0(g 



2 7 (1) 



Eq. < J3 -3[ ) means in terms of these functions 

Iia(6,xo/L) = 0(a/L) , V0,xq/L, at amo = am c 



(5.3) 



(5.4) 



(5.5) 



Employing this condition at tree level one finds that the critical mass m c vanishes, i.e. am c = 0(g 2 ). 
Here we determine the first order coefficient in 



am c am>^ g 2 + 0(g 4 ) 



Using expansions ( p.l| ) and ( p.q ) in ( p.3| ) gives 

h A = ho + g 2 [h x + am^ h 2 + (z^ + 2z[ 1} ) ho) + 0(g 4 

with 



h = Ldof^ (xq) , hi = Ldoffl (x ) , h 2 



damo 



ho, 



(5.6) 



(5.7) 



(5.8) 



all defined at amo = 0. The tree level amplitude ho turns out to be 0(a 2 ), whereas h 2 diverges 



linearly in a. Therefore ( p.5| ) amounts to 



am 



(i) 



-hi/h 2 + 0(a 2 ). 



(5.9) 



Since the amplitude h\ 9 as fj^\ is linear in the coefficients c\ introduced in ( pT2| ) the finite and 
the 0(d) part of —h\jh 2 yield two linear equations for each pair of 9 and xq/L. But the limit 
a — > 0, defining the critical mass, turns out to be independent of and xq/L giving a single relation 
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between the critical mass and the coefficients. The 0(a) part has to vanish, thus demanding cs = cp. 
In terms of the coupling constants the result is 

am c = -0.7698004(1) x ((2N - l)/2g 2 + g 2 /2-g 2 ) +0(g 4 ) , (5.10) 
g 2 P /g 2 s = l + 0(g 2 ). (5.11) 

To compare with the large N result of [Q] we rescale gf — > gj/N and take the N — > °o limit in ( 5.10|) 



am c = -0.7698004(1) x g| . (5.12) 

This is the full large N result. The authors neglect the vector 4FI from the beginning, but their 
result is not changed if it is taken into account. In the large N limit gp is a nontrivial function of g 2 s 
which is reproduced by our result at first order. 

6. Conclusion and Outlook 

Self-coupled fermions in two dimensions with U(N) flavour and continuous chiral U(l) sym- 
metry are asymptotically free in the large N limit and in perturbation theory. We establish this the- 
ory on the lattice with the Wilson discretization. In general there are three dimensionless coupling 
constants and the bare mass parameter due to the explicitly broken chiral symmetry. Continuum 
Ward-Takahashi identities are used to enforce chiral symmetry in the continuum limit of the lattice 
theory fixing the parameters of the action. At first order in perturbation theory, when all couplings 
are small, we determine the critical mass and the ratio of the pseudo scalar and the scalar coupling. 
The vector coupling is not fixed at this order. It is not needed for the renormalisation of the ratio. It 
appears in the expansion of the critical mass which is not seen at large N. Nevertheless our result 
contains the full large N critical mass and also the ratio of the couplings is correctly reproduced. 
A next to leading order calculation is almost finished and covered in a forthcoming publication 
With this result we will be in the position to calculate a renormalised coupling and thus a step 
scaling function in the controlled environment of Wilson fermions. Then a direct comparison with 
other fermion actions like Ginsparg-Wilson and staggered is possible and will be done. 

Acknowledgment We would like to thank Rainer Sommer and Francesco Knechtli for helpful 
discussions and Ulli Wolff for reading the manuscript. 
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